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1 INTRODUCTION 


1 . 1 Statement of problem 

It is required to somehow suspend or levitate an aerodynamic 
model in such a way that a large clear air volume can be 
maintained around the model. It has been decided that this could 
be achieved most satisfactorily if all levitation/suspension 
hardware were positioned on or behind a single plane surface, at 
some distance from the center of the model. This distance would 
be of the order of 3 to 4 feet. Further, there appears to be some 
advantage in choosing the floor to be the plane surface, forming 
a Magnetic Levitation (ML) system (Fig. 1.1). This report 
represents a preliminary assessment of the technical feasibility 
of such a system. 

The model's own weight will be the principal force acting on 
the . levitated model. Rather small aerodynamic or other loads are 
anticipated. The model's spatial position and orientation are to 
be maintained to high accuracy, with very small allowable 
undemanded motion around any particular location. The position 
and attitude must, however, be variable, under operator control, 
preferably over a wide range. The rate of movement from one 
position/attitude to another may be relatively slow. 


1.2 Background - wind tunnel MSBSs 

Magnetic Suspension and Balance Systems for wind tunnel 
models (MSBS) have been under development for nearly 30 years. 
Around 15 systems have been constructed by ten different 
institutions, worldwide (Ref.l). Only four systems are presently 
known to be operational : 

13 inch NASA Langley Research Center (formerly AEDC) 

6 inch MIT/NASA Langley Research Center 

7 inch University of Southampton 

3 inch Oxford University 

General features common to all known systems include small 
physical scale, the 13" LaRC system being the largest yet 
constructed, and relatively low magnetic force capability, 
typically less than 5 times the weight of the model ' s magnetic 
core. With the exception of the second University of Virginia 
system, all have utilized conventional copper electromagnets, 
with or without forced cooling. Two systems (Southampton and 13" 
LaRC) are equipped with digital control systems, otherwise simple 
analogue controllers have been used. A wide variety of techniques 
for sensing the position and attitude of the suspended model have 
been employed, including many types of optical system, 
electromagnetic sensors and X-rays. Still more techniques have 
been proposed, though no system currently in operation can be 
regarded as entirely satisfactory. Model magnetic cores have been 
soft iron or permanent magnet, though a superconducting solenoid 
has been demonstrated. 
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Design studies of Large MSBSs (LMSBS) have been performed 
recently, (Refs. 2,3,4) with the important conclusion that 
systems for wind tunnels up to 8 foot test section size are 
feasible with essentially existing technology, though potentially 
costly to construct. 


1.3 Comparative technical difficulty of MSBS and ML 

A comparison of this type is appropriate and illuminating. 

Physical scale 

Although existing MSBSs are all relatively small, the design 
studies presently available give a high degree of confidence that 
large scale facilities (up to 8-foot test section size) could be 
successfully constructed. The single major technical change 
required for the larger scale systems, compared to existing 
MSBSs, would be the replacement of the usual copper 
electromagnets with superconducting electromagnets. This is for 
two important reasons: 

1) Creation of the relatively intense magnetic fields 

required to oppose the high aerodynamic loads on the 
model . 

2) Power economy. The steady state power consumption of 

copper electromagnets increases rather rapidly as 
system size increases and can become outrageous. This 
problem is perhaps the more relevant for ML. 

Force and moment capability 

MSBSs have to be capable of magnetically opposing 
uncontrolled and unknown aerodynamic loads on the model. These 
loads can be large with respect to the model's mass and will be 
so in LMSBSs. The absence of these loads in the ML application 
has two important effects : 

1) Reduction of required intensity of magnetic fields and 

field gradients. This apparent advantage is, however, 
offset to some degree by the restrictions on location 
of the levitation electromagnets. 

2) Considerably reduced dynamic force and moment 

requirements. MSBSs have to be able to counter rapidly 
changing aerodynamic loads and usually are capable, as 
a result, of driving the model through operator 
demanded, oscillatory motions of significant amplitude 
and frequency. Quasi-steady, undisturbed levitation is 
possible with substantially reduced dynamic capability, 
notably in the area of electromagnet power supplies. 

Position/attitude precision 

The required precision for MSBS and ML appear to be somewhat 
similar, specifically of the order of thousandths of an inch in 
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position and hundredths of a degree in orientation. This has been 
achievable for many years with MSBSs, albeit under controlled 
conditions. Long-term drift has been the most difficult problem, 
almost entirely related to position sensing techniques. New 
approaches in this area promise to eliminate this difficulty. 

A model in an MSBS can never be absolutely stationary, since 
the systems are open-loop unstable and rely on position and 
attitude error feedback for stabilization. The exact amplitude of 
motion of the model around a set-point is a strong function of 
the amplitude of the disturbances acting on the model 
(aerodynamic etc.), the quality of electromagnet power supplies, 
the resolution of position sensors and the sophistication and 
adjustment of the control system. With wind-off (low 
disturbances) , extremely steady suspension has been observed at 
the University of Southampton (certainly +- 0.001" position). 
Fairly good results have been measured with the 13" LaRC MSBS (+- 
0.003" position) despite several recognized shortcomings in 
hardware ( Fig . 1 . 2 ) . 


Position/attitude range 

MSBSs have mostly been restricted to a rather small range of 
model positions and attitudes. The exact ranges have generally 
been set by position sensor capability, but there are fundamental 
restrictions caused by the gross variation in magnetic couplings 
between electromagnets and model that occur with large angular or 
translational motion. This difficulty has been addressed at some 
length in Ref. 5 and limited but successful demonstrations have 
been made of extreme attitude capability (Fig. 1.3). Extensive 
further work in this area is required and some is understood to 
be underway at Southampton. It is, however, confidently stated 
that there is no fundamental limit to model position or attitude 
capability, provided that : 

1) The correct magnetic field and field gradient components 

can somehow be created in the appropriate strengths and 
combinations . 

2) The extensive couplings between model degrees of freedom 

that occur due to the powerful and variable magnetic 
couplings can somehow be accomodated. 

3) The position, attitude and trajectory of the model can be 

accurately monitored. 

Control systems 

MSBSs utilize multi-loop control systems of the general form 
shown in Fig. 1.4. Control algorithms may be implemented 
digitally or by analogue circuitry. This general form is 
certainly applicable to ML. The compensation algorithms (series 
lead-lag) universally used to date are somewhat antique, but can 
be made to perform well with some care in adjusting gain and 
frequency parameters . It is felt that ML can be expected to 
operate successfully with this form of algorithm and, with a 
digital control system, new algorithms can easily be installed 
whenever they become available. 
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The "prefilter" and "translator" functions will not be 
trivial for ML. They exist so that the controller can operate in 
model degrees of freedom, strongly felt to be the best approach. 
If the model is to be capable of large attitude or position 
excursions, then the prefilter and translator become large 
matrix operators with variable coefficients. The general theme 
explored in Ref. 5 was that these coefficients could be pre- 
calculated to adequate accuracy and adjusted in real-time by 
using look-up tables or programmed functions or a combination 
thereof in a control computer. This has never been demonstrated 
experimentally, though suspension at various discrete attitudes 
was acheived with pre-calculated de-coupling/coupling terms in 
the prefilter and translator. 

Data acquisition 

In the wind tunnel application, whole-body aerodynamic force 
and moment data may be acquired by monitoring the electromagnet 
currents and comparing to previous calibrations. Lack of 
efficient calibration techniques has been a long standing and 
unresolved difficulty. This kind of calibration and data 
acquisition is not required in the ML application. 

Electromagnet power supplies 

Power supplies have been a source of many problems in MSBS 
development (Power supply is actually something of a misnomer, 
the proper description being power amplifier, since controlled 
currents must be supplied to the load electromagnet by varying 
the amplifier output voltage) . Current and voltage requirements 
can be rather high, many kVA even at small scale, and frequency 
response must extend to a few tens of cycles per second. The load 
current must be relatively smooth for high quality control of 
model position although the high inductances typical of large 
electromagnets is of assistance here. These are demanding 
specifications but can be met by power supplies intended for 
large D.C. servomotor drives. Units are known to be commercially 
available up to around 500 kVA rating (1000 Amps, 500 Volts) . The 
most difficult requirement for MSBSs has been for bipolar 
operation. This arises partly since the weight of the model can 
be relatively small compared to other loads, thus the 
electromagnet currents required merely to suspend the model 
(wind-off) are low. Further, as the model is maneuvered over a 
range of attitudes, large changes, including sign reversals, of 
electromagnet currents are encountered. Since the model's weight 
is the dominant force in ML, monopolar current capability, at 
least in some electromagnets, may prove to be adequate, though 
the magnetic coupling changes with attitude and position 
variations require careful study in this regard. 

Magnetic configuration 

Advanced MSBS designs tend to develop as a tightly packed 
array of many electromagnets surrounding the wind tunnel test 
section. The geometry and configuration of these electromagnets 
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is crucially important for several reasons : 

1) Creation of the required field and field gradient 

components . 

2) Achievement of 1) with high efficiency 

3) Achievement of 1) with predictable and manageable 

magnetic couplings over the required range of model 
positions and attitudes. 

The same requirements apply to ML systems, though the 
designer's freedom to locate electromagnets around the model is 
seriously reduced. 


Summary 

With the sole exception of magnetic configuration, it is 
seen that ML is of comparable or lower technical difficulty than 
LMSBS. Synthesis of an appropriate magnetic configuration is 
essentially an intellectual problem, to be tackled and resolved 
in the early phases of design. The technical risk for an ML 
system, meeting the broad specifications and operational 
requirements presently considered, is therefore judged to be 
relatively low. 
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Fig. 1.1 Demonstration magnetic levitator - conceptual 
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Fig. 1.2 Demonstrated steadiness of model suspension 
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Fig. 1.3 High angle-of-attack suspension with MSBS 



Fig. 1.4 Typical MSBS control loop 
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2 GENERAL DESIGN APPROACH 


2 . 1 Choice of model core 


Three options are available for the magnetic core of the 
levitated model: 

1 ) Permanent magnet 

2) Soft-iron 

3) Persistent superconducting solenoid 

All have been used to various degrees with MSBSs. Soft-iron 
cores have been utilized in two fundamentally different ways, one 
being operation in an attraction mode (eg. 13" LaRC) , where the 
required magnetization of the iron is naturally achieved by the 
applied fields, as shown in Fig. 2.1. Alternatively, the 
magnetization can be artificially introduced by quasi-uniform 
applied fields (eg. 6" MIT/LaRC) , also illustrated in Fig. 2.1. 

This is normally only successful if the subsequently applied 
fields and field gradients exhibit a certain degree of symmetry, 
otherwise the induced magnetization may be lost. 

The requirement to levitate the model above a plane floor in 
this application suggests operation in a repulsion mode. Further, 
there appears little possibility of elegantly achieving the 
required uniformity of magnetizing and symmetry of other fields 
from the restricted available locations of electromagnets. Option 
2) above is therefore dismissed. 

For large scale applications, a persistent superconducting 
solenoid has some advantages. The high levels of effective 
magnetization that are attainable would result in significant 
economies in the sizing and powering of the levitation 
electromagnets. A solenoid of representative scale has already 
been successfully operated in a small MSBS (Ref. 6). Technical 
difficulties in the use of superconducting solenoids in MSBSs 
include limited run-time (of the order of 1 hour) , helium 
venting, calibration uncertainties and mechanical design to 
sustain high loads. The levitation application does not require 
high load capability, nor force and moment calibration and the 
residual difficulties of run-time and helium venting are regarded 
as relatively minor development problems. Unfortunately, the 
performance of superconducting solenoids is very sensitive to 
scale, improving rapidly, relative to alternative core materials, 
with increasing scale. 

For the scale of device presently envisaged, the appropriate 
choice of core material is felt to be permanent magnet. Powerful, 
stable materials are available (Rare-Earth Cobalts) in large 
sizes (at least 2" dia. by 0.5" length, stackable). The control 
and dynamics of a superconducting solenoid are virtually identical 
to a permanent magnet model, thus a superconducting solenoid 
could be utilized later, providing a performance upgrade of some 
sort, if required. 

At the present time, representative properties for Rare- 
Earth Cobalt materials are taken to be: 
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Magnetization (Saturation Induction) 
Density 

Coercive force (greater than) 

Magnetic stability 
Mechanical properties 


1 Tesla 
8000 kg/m 3 
800 kA/m 
Very good 
Hard, brittle 


The only problems that are likely to be encountered are the 
mechanical brittleness and the the fact that the magnet cannot be 
"turned off", presenting certain cleanliness and housekeeping 
problems. 


2.2 Magnetic configuration 

The general requirement has been concluded as being that of 
levitation of a model in "normal" orientation (that is, major 
axis horizontal) by an array of electromagnets buried beneath a 
plane floor (Fig. 1.1). This configuration is expected to be 
statically unstable (that is with steady magnetic fields) in 
several degrees of freedom. Vertical translation and pitching 
rotation may be stable, provided other motions are suppresed. 

Roll rotation may be neutrally stable with axisymmetric models or 
positively stable with the model center of gravity set lower than 
the magnetic center. Now, some open-loop unstable degrees of 
freedom are always present in magnetic suspension systems of this 
general class. Further, any stable degrees of freedom tend to 
exhibit unsatisfactorily low damping. Feedback control of all 
degrees of freedom, providing artificial stiffness and damping, 
is the normal and correct approach, thus the extensive open-loop 
instabilities of this configuration are of no particular concern. 

It is desired to control the model over a significant 
range of attitude, perhaps : 

Pitch : +- 20 degrees 

Roll : +- 20 degrees 

Yaw (Azimuth) : 360 degrees 

Since the levitated model is part of a control loop, some 
angular motion is naturally available. Experience with MSBSs 
indicates that the pitch and roll requirements shown above might 
be met with little difficulty. The azimuth requirement is 
significantly more challenging. Presuming that the model position 
and attitude can be accurately monitored over the required range, 
itself not a trivial task, the synthesis and verification of 
appropriate electromagnet configurations and development of 
decoupled control algorithms still remain as serious 
difficulties. These points are discussed more fully in Section 5. 


2 . 3 Other hardware 


Sufficient experience already exists in the use of digital 
control systems with MSBSs to conclude that a digital controller 
is the correct choice for ML. The control task can be considered 
as being divided into three main sections (after Ref. 7): 
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1) Stabilization 

2) Operator interface and supervisory functions 

3) Decoupling of position sensor data and power supply 

demands 

The computational capacity required for 1) above can be 
reliably extrapolated from current MSBSs, where PDP-11 class 
minicomputers are successfully used. The operator interface and 
supervisory functions are not computationally intensive, but have 
proved somewhat awkward to integrate into control software. A 
separate processor might be the best choice here. At the present 
time it is difficult to estimate the computational capacity 
required for the decoupling task. It is thought likely, however, 
that equivalent or even greater capacity might be required for 
this task as is necessary for stabilization. Since these two 
tasks would have to be fully integrated, a small machine of the 
VAX class might be the minimum appropriate choice for the main 
control computer. 

Solid-state photodiode array cameras are under investigation 
as model position sensing devices for application to LMSBSs 
(Refs. 2, 8) and appear to be an excellent choice for ML. The 
advantages of these devices, compared to known alternative 
techniques, are ideally matched to ML requirements : 

1) High achievable accuracy 

2) Near-absolute repeatability 

3) Insensitivity to magnetic fields 

4) Wide posit ion/ attitude range 

The potential disadvantages for LMSBS application are either 
irrelevant here or soluble : 

1) Light refraction with variations of air density 

2) Significant computational overhead (image processing) 

3 ) High cost 

Some commercial ^position sensing devices of this general 
class are available . It seems certain that a system meeting ML 
requirements can be developed. Choice of optical "target" on the 
levitated model and provision of high-speed real-time processing 
hardware are seen as the most serious difficulties. Integration 
of the image processing tasks into the main control computer is 
strongly discouraged. 


* 


Saab-Scania, 


Eloptricon Division, 


"Eloptopos" 
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a) Model magnetization achieved by suspension field 



b) Model magnetization achieved by quasi-uniform applied field 


Fig. 2.1 Model magnetization options 
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3 MAGNETIC THEORY 


3 . 1 Force and moment production 


This area has been covered at some length in Refs. 5 ,9, 10 and 
elsewhere and will only be reviewed here in terms specific to the 
ML application. 

The governing equations for the force and moment on a 
magnetically suspended/ supported model are as follows : 


F = yU 0 J M.VH dV T + rX(M.VH) dV 


-where M is in Amperes/metre. If a permanent magnet model 
core is used then M is known and constant. H will be the 
"external" H created by the electromagnet array. For preliminary 
design, particularly where the model core is relatively small 
with respect to the suspension electromagnets, the integrals may 
be eliminated by assuming that the relevant field or field 
gradient components are uniform over the volume of the core : 


F = ^ 0 V(M.VH 0 ) T = yU 0 V(MXH 0 ) + second order terms 

Hq is evaluated at the model core centroid. >CLM is frequently 
replaced by J (polarization, Tesla) if the core is a permanent 
magnet . 


Defining axes in the usual aircraft sense (Fig. 3.1) and 
assuming a conventional axially magnetized core, it is 
straightforward to estimate the magnetic field required to 
support the core's own weight. 

—magnetic “ ’^gravity = - pgv k = J x VH Jtz k 

Using representative permanent magnet material data from Section 
2.1, the required field gradient is found to be : 

H _ = 78,480 A/m/rn = B = 0.0986 T/m 

xz xz 

Since there is some deadweight to be carried by the core and 
since an excess force and moment capability is required for 
control and maneuvering, a Load Factor (n) , defined in the normal 
aircraft sense is introduced : 


—magnetic n * ^ "^gravity ^ 

At this point it is assumed that a reasonable design value for 
Load Factor is approximately 2, implying that the core could 
support an aerodynamic envelope perhaps 50% of its own weight, 
still leaving around 30% overcapacity for control and 
maneuvering. Refined values for Load Factor could be determined 
from comprehensive system simulation studies. Extra capability 
that may be required in large position or attitude excursions 
must be dealt with separately. 

Criteria for axial (x-wise) and lateral (y) force 
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capabilities are similarly arbitrary at this stage, but have been 
chosen to be equivalent to a Load Factor of 0.2, or 10% of the 
vertical (z) force capability. Thus : 

H xx “ H xy “ 15,700 A/m/m = B xx = B xy =0.02 T/m 

Torque or moment production has not generally been a design 
driver for MSBSs apart from the roll degree of freedom, which 
will be dealt with separately. For an assessment of the order of 
magnitude of magnetic field required, consider that the model 
core's mass and magnetic centers may be offset from each other by 
some amount. The force due to gravity and the opposing magnetic 
force, for instance, will now create a torque of some value. If 
typical model core dimensions are taken to be 9" long (22.5 cm) 
and 2" diameter (5 cm) and if the offset of magnetic and mass 
centers was 20% of the length of the core, then the maximum 
torque would be : 

-max ■ 2 f* gV *° • 225*0.2 Nm = J x VH z 2 

Solving for H gives a value of 7063 A/m/m, or B of 0.009 Tesla. 
A value for H z can be derived similarly. z 

Generation of rolling moment has traditionally been a 
difficult problem with MSBSs. Refs. 5 and 11 contain extensive 
discussion of a variety of "D.C." techniques, that is using 
steady applied fields and non-uniform or non-axisymmetric model 
core magnetization distributions. For the time being, it can be 
assumed that some magnetic roll torque will be required, though a 
levitation system might be operable with very low roll torque 
capability, since it is possible to arrange for that degree of 
freedom to be naturally stable (a frequent practice with MSBSs) . 
Positive stiffness is achieved by positioning the model's mass 
center below its magnetic center, as shown in Fig. 3.2, either by 
appropriate magnetic design or by ballasting. Magnetic roll 
torque is then only required for the introduction of artificial 
damping and for control of the roll datum orientation. 


3.2 Magnetic couplings 

As previously discussed, extreme position or attitude ranges 
require that the problem of the variation of magnetic couplings 
be addressed. The theoretical background relevant to attitude 
variations is summarized here, from Ref. 5. Coupling variations 
with model translations have not been a great concern for MSBSs, 
where the model remains roughly centered in the test section for 
all normal testing, and has not been studied at any length. It is 
felt, however, that this aspect of the problem is no more 
difficult than that arising during attitude changes. Indeed, 
coupling variations arising during translations where the 
model/levitation system retains extensive symmetry are likely 
only to require control loop gain scheduling. 

Using conventional Euler angles for specification of model 
orientation (Fig. 3.1), it can be shown that field components in 
"model" axes are related to those in "levitation" axes as follows : 
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Thus the simplified torque equation becomes : 

T 7 ® VM ' X (AH) + second order terms 

Field gradient components in model axes relate to those in 
levitation axes thus : 
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The simplified force equation becomes : 

F* = V(M'.AV)AH 

It is seen that field components behave as vectors during axis 
rotations, whereas field gradient components do not, and that the 
couplings from applied field gradient components (levitation 
axes) to forces in model axes will be extremely complex. 

If individual or clearly identifiable groups of 
electromagnets created, or could be arranged to create, single 
field gradient components, then the couplings described above 
could be accomodated straightforwardly by including an "inverse" 
coupling matrix of some sort in the control system. The important 
and crucial difficulty arises due to the fact that individual 
electromagnets create combinations of several, or all, field 
gradient components simultaneously, except under rather specific 
geometric conditions (eg. Helmholtz pair) . Further, there are 
usually more electromagnets available than there are model 
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degrees of freedom to be controlled. Under these conditions there 
will not be a unique solution for decoupling. Identification 
of the optimum solution (yielding greatest force and moment 
capability for a given set of electromagnets) presently involves 
repetitive calculation (Ref. 5). 

It was presumed in Ref. 5 that decoupling matrices could be 
satisfactorily precalculated, using computed or measured magnetic 
field distributions. Stable suspension was achieved with the 
Southampton University MSBS in this manner, though only using 
fixed decoupling coefficients. In the same way that some position 
and attitude capability is available about the conventional datum 
position and orientation in MSBSs, some position and attitude 
capability exists around each design decoupled attitude. The next 
logical step would be to precalculate couplings over a wide range 
of positions and attitudes and to arrange for the decoupling 
matrix to smoothly transition between known forms as the 
levitated model is maneuvered. It is felt that this approach will 
be successful, but experimental demonstration and verification is 
certainly required. 


15 








4 PRELIMINARY MAGNETIC DESIGN and OPTIMIZATION 


4.1 Simple , uncontrolled levitation 

One of the simplest conceivable levitation arrangements is 
that of two circular coils as shown in Fig. 4.1. It is presumed at 
this point that if these coils were arranged to take the weight 
of the model then several (smaller) coils could be distributed in 
some fashion so as to provide the necessary force and moment 
capability in other degrees of freedom. This would represent a 
baseline design for an ML system. The required sizes and 
specifications of the levitation coils are therefore of interest. 
Using the computer program FORCE (Appendix 1) , it is 
straightforward, though somewhat tedious, to compute the 
performance of families of coil pairs of different sizes, 
covering a range of possible geometries (Fig. 4.1). This procedure 
requires that the coil current density be specified. 


4.1.1 Choice of superconductor specifications 

A value of 1500 Amps/cm 2 is chosen for coil current density. 
This is consistent with previous Large MSBS designs (Refs. 2, 3, 4) 
and is felt to be reasonable and conservative (cryostable 
conductor) . The figure represents a mean current density across 
the whole coil cross section, including some structure and inter- 
winding insulation, but not including the helium dewar or any 
external structure. Adequate inter-coil clearances must be 
allowed for these components. A minimum spacing of 5 centimetres 
between adjacent coil windings has been assumed, based loosely on 
existing LMSBS design data. 

Higher current densities are possible with small coils 
(adiabatically stable conductor) , up to ten times the figure 
quoted above. Assessment of the feasibility of this type of 
conductor in this application would, however, require expert 
opinion and will not be addressed further here. 

Calculation of coil inductances is necessary for sizing of 
the power supplies, but requires that the conductor cable current 
be specified. Superconducting cables are generally designed for 
specific current levels and development of a new conductor is 
likely to be a lengthy and risky process. Large MSBS designs have 
therefore concentrated on the use of off-the-shelf conductor 
designs, specifically an llkA cable developed for Argonne 
National Laboratory (Refs. 3, 4) and a 50kA cable developed for Los 
Alamos (Ref. 2). Both of these appear to be too high a current 
rating for the present application, therefore an arbitrary cable 
current of 1500 Amps is used extensively in this report (1 cable 
per square centimetre) . Adjustment of calculated inductances to 
suit specific conductor specifications is straightforward. 


4.1.2 Two-coil levitation system 

Computed B field gradient components generated at the 
xz 
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model center by pairs of electromagnets of specified geometry are 
summarized in Fig. 4. 2 and the Tables below. 


• 1 
Q I 
• 1 
O 

(m) I.D.(m) X (m) 

Y (m) 

B xz (T/m) 

Vol. (m 3 ,each) 

Square cross-section windinqs (Winding A) 


1.0 

0.5 

0.525 

0.25 

0.0598 

0.1473 

1.4 

0.7 

0.725 

0.35 

0.1734 

0.4041 

1.6 

0.8 

0.825 

0.4 

0.2509 

0.6032 

1.8 

0.9 

0.925 

0.45 

0.3387 

0.8588 

2:1 

Flattened 

coil (Winding B) 



1.4 

0.7 

0.725 

0.175 

0.1044 

0.2020 

1.6 

0.8 

0.825 

0.2 

0.1522 

0.3016 

1.8 

0.9 

0.925 

0.225 

0.2069 

0.4294 

2.0 

1.0 

1.025 

0.25 

0.2667 

0.5890 

4:1 

Flattened 

coil (Winding C) 



1.6 

0.8 

0.825 

0.1 

0.0845 

0.1508 

2.0 

1.0 

1.025 

0.125 

0.1490 

0.2945 

2.4 

1.2 

1.225 

0.15 

0.2227 

0.5089 

2.8 

1.4 

1.425 

0.175 

0.3003 

0.8082 

1:2 

Thickened 

coil (Winding D) 



1.0 

0.5 

0.525 

0.5 

0.0900 

0.2945 

1.2 

0.6 

0.625 

0.6 

0.1606 

0.5089 

1.4 

0.7 

0.725 

0.7 

0.2514 

0.8082 

1.6 

0.8 

0.825 

0.8 

0.3591 

1.2064 

Electromagnet 

inductances 

can be 

estimated using 

formulae in 


Ref. 12 and emerge as follows : 


L, - 1.274 * Outside radius * N 2 MicroHenries 
A 

L- - 1.523 * Outside radius * N 2 MicroHenries 
B 

L n = 1.687 * Outside radius * N 2 MicroHenries 

v 

L d - 0.959 * Outside radius * N 2 MicroHenries 


Only one size of each geometry of winding produces the required 
field gradient (see Fig. 4. 2). Identifying these sizes and 
assuming a 1500A cable current, the following data can be 
calculated (all data per coil, see also Fig. 4.3) : 
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Winding O.D.(m) Vol.(m 3 ) No. turns Induct. (H) Energy (MJ) 


A 

B 

C 

D 


1.47 0.468 

1.78 0.415 

2.28 0.436 

1.29 0.632 


1351 

1.708 

990 

1.179 

812 

1.269 

2080 

2.675 


1.921 

1.326 

1.428 

3.009 


The electromagnet sizes calculated above seem to be well 
within the bounds of present technology. Numerous electromagnets 
of larger sizes, volumes, inductances and stored energies have 
been constructed for a variety of applications. 

Power supply requirements are a direct function of 
inductance (or stored energy) . Significant dynamic current 
capability is necessary for control of the levitated model. The 
required dynamic current capability is not, however, clear since 
it is also a function of the level of disturbances encountered 
and the characteristics of the control loop, particularly 
suspension stiffness. Based on experience with MSBSs, a "ball- 
park" figure can be established. Again, refined estimates could 
be obtained from extensive system simulation studies. If a 
required variation of current is chosen to be plus-or-minus 10% 
of maximum available current, at a frequency of 1 cycle per 
second, then the required voltage for near-optimum electromagnet 
geometry is given by: 


v max * L “ L ■ 1.2*300/0.5 - 720 Volts 


Power supplies of this order of performance are commercially 
available. Individual electromagnets can be subdivided and fed 
from a number of smaller power supplies if required, but there 
are some disadvantages in this approach. Low voltage, D.C. 
supplies for sections of electromagnets have been used with MSBSs 
("Bias" supplies) but are regarded as an unattractive option 
here, due to the inductive couplings between adjacent windings. 

Optimization criteria are to seek a minimization of 
winding volume in order to economize on the use of superconductor 
together with a minimization of inductance (or stored energy) so 
as to minimize power supply requirements. Fig. 4. 3 clearly shows 
that this is achieved for this phase of design with relatively 
flat electromagnets. The variation of important parameters is 
quite weak around the optimum design, this being the usual case. 

A winding aspect ratio of around 2 appears to be close to 
the local optimum. Specifications for electromagnets of this 
geometry would be : 


Outside diameter 

Inductance 

Stored energy 

Turns (of 1500A cable) 

Length of conductor 

Mass < 


1.78 m 
1.18 H 
1.33 MJ 
990 

4150 m 
3000kg 


Field at model centroid B xz 


0.2 T/m; B x - 0.113 T 
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4.1.3 Copper electromagnet feasibility 

At this point, the feasibility of a conventional copper 
electromagnet system can be addressed. The author is not aware of 
any "magic number" for usable current density for copper 
windings. Quite high current densities are achievable with 
careful design and effective cooling techniques. Indeed, the 
technical requirements for cryostable superconducting 
electromagnets, where significant lengths of normally conducting 
cable must be re-quenched by the liquid helium coolant, indicate 
that current densities of the same order as 1500 A/cm 2 could be 
achieved. If 75% of the cross section of the windings were copper 
(remainder for turn-to-turn insulation and coolant passages) then 
it is easy to calculate the steady state power consumption for 
the optimum winding geometry identified above. 


Outside diameter 
Inside diameter 
Depth 

Current density 
Packing factor 
Power at max. current 


1.78 m 

0.89 m 

0.445 m 2 

1500 A/cm 2 overall 

0.75 

2.12 MW 


The latter figure represents the steady-state electrical power 
demand, some fraction of required power supply MVA and, depending 
on cooling technique, the refrigerating power to be delivered to 
the load. It is felt that this power consumption would be 
unreasonable for the type of facility presently contemplated. 

It must be emphasized, however, that a copper system is 
probably technically feasible. Detail optimisation of the copper 
electromagnets is undoubtedly possible, since reduction of 
current density is desirable for reduction of power dissipation. 
No further analysis will be attempted here. 

Iron cores are frequently used with MSBS electromagnets in 
order to enhance the field generated by any given size of 
electromagnet. However, iron cores are not used extensively with 
superconducting electromagnets since the relatively intense 
fields in the bore of the windings tend to completely saturate 
the iron, reducing its effectiveness. Further, the iron must 
either be within the liquid helium dewar, where it would 
represent an extra thermal mass to be cooled, or be external, 
thus separated from the windings by structure and dewar walls. 
Calculation of the effectiveness of iron cores is extremely 
difficult (see Appendix 1) and will not be attempted here. Since 
the two levitation electromagnets studied so far would be 
essentially monopolar devices and of opposite polarity relative 
to the model location, the use of some form of iron yoke as a 
flux return path (Fig. 4. 4) may be an attractive option and should 
be studied further. 
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4.2 Controlled, single point levitation 


4.2.1 W 2+6 W electromagnet configuration 

As previously mentioned, the two-electromagnet levitation 
system can, in principle, be extended to a system with full 
control capability in all degrees of freedom by the addition of a 
number of smaller electromagnets. A relatively simple arrangement 
of six additional electromagnets is shown in Fig. 4. 5. These 
electromagnets would be required to generate at least the 
following fields or field gradient components : 

B (Yaw moment) 

y B (Sideforce) 

** B xx (Axial force) 

Rolling moment will be dealt with separately. B (Pitching 
moment) can be generated by the two levitation electromagnets , 
though the additional control electromagnets might be required 
for trimming. The effectiveness of various sizes of additional 
electromagnets has been calculated : 


Geometry 

Y (m) 

O.D. (m) 

I.D. (m) 

Depth (m) 

3 

Volume (m ) 

AA 

0.8 

0.55 

0.275 

0.1375 

0.0245 

BB 

0.9 

0.686 

0.343 

0.1715 

0.0475 

CC 

1.0 

0.83 

0.415 

0.2075 

0.0842 

DD 

1.1 

0.982 

0.491 

0.2455 

0.1394 

EE 

1.2 

1.138 

0.569 

0.2845 

0.2170 

FF 

1.3 

1.3 

0.65 

0.325 

0.3235 


Current 

sense 

B„ v 



Geometry 

1 2 

3 4 

xy 



AA 

+ 

+ 

0.001 



BB 



0.0026 



CC 



0.0052 



DD 



0.01 



EE 



0.0172 



FF 



0.0276 



Geometry 

Current 
1 2 

sense 
3 4 

B z 

B xx 

B B 

yy zj 

EE 

+ + 

+ + 

-0.0022 

0.0158 

0.0008 -0.1 

FF 



-0.0038 

0.0226 

0.0034 -0.1 


It is seen that the required components B xy and B xx can be 

generated by electromagnets of geometry comparable to EE and FF 
above. The "stray" field components accompanying the B xx 

component are of no great concern, since B^y and B zz have no 
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primary effect with the core geometry chosen. B z would cause a 

pitching moment, but this is easily countered by the levitation 
electromagnets, or electromagnets 5 and 6 from Fig. 4. 5. 


Current sense B 

Geometry 5 6 


B B 

xx yy 


B 


zz 


EE + + -0.0178 0.0286 -0.0296 0.0008 

FF -0.0238 0.0398 -0.0462 0.0066 

Considerable scope for optimization exists, since there is a 
clear mechanism for trade-off of winding volume between the two 
levitation electromagnets and the six control electromagnets. As 
the diameters of the levitation electromagnets are reduced, so 
the control electromagnets appear to be able to move into more 
favourable locations. Circular electromagnets are preferred by 
manufacturers, due to mechanical simplicity and relatively 
uniform winding stresses, though study of elliptical or racetrack 
electromagnets might be worthwhile. 


4.2.2 "e+O" electromagnet configuration 

Since the six electromagnets added for control and 
stabilization in Section 4.2.1 became of comparable size to the 
two levitation electromagnets, it may be reasonable to eliminate 
these two and bring the remainder into a more tightly packed 
cluster, enlarging them so as to satisfy the levitation 
requirement. The basic configuration is shown in Fig. 4. 6 with 
calculated levitation performance shown in the tables below. 

Examination of this data shows that this configuration is 
relatively unpromising. Even rather large electromagnets fail to 
provide sufficient field gradient for levitation (0.2 Tesla/m). 

The first obvious approach in an optimization procedure 
would be to reduce the diameter of the two control electromagnets 
(2,5 in Fig. 4. 6) so as to allow the other four to be moved into 
more favourable locations. 


O.D. (m) 

I.D. (m) 

Depth (m) 

3 

Volume (m ) 

B X (T) 

B xz ( T / m ) 

Square cross-section windings 




1.2 

0.6 

0.3 

0.254 

0.0626 

0.0706 

1.6 

0.8 

0.4 

0.603 

0.1168 

0.083 

2.0 

1.0 

0.5 

1.178 

0.1726 

0.0738 

2.4 

1.2 

0.6 

2.036 

0.225 

0.052 

Thickened 

coil 





1.0 

0.5 

0.5 

0.294 

0.0668 

0.0954 

1.2 

0.6 

0.6 

0.509 

0.1062 

0.123 

1.4 

0.7 

0.7 

0.808 

0.1522 

0.1444 

1.6 

0.8 

0.8 

1.206 

0.2016 

0.1566 

1.8 

0.9 

0.9 

1.718 

0.304 

0.1588 
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Field 


x 
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Inductance (Henries) 
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Fig. 4. 5 tl 2+6 lt electromagnet configuration 



Fig.4.6 ”6+0 11 electromagnet configuration 




5 FURTHER MAGNETIC DESIGN 


5.1 Roll control 


In model axes, several field or field gradient components 
are already primary producers of force or moment components : 

B Yaw moment 

y B Pitch moment 

Z B w Axial force 

xx B Sideforce 

™ B xz Vertical force 

An effective roll moment technique would utilize some other 

field or field gradient component so as to remain uncoupled, to a 

first order, from other degrees of freedom. The axial field 

component (B ) cannot produce roll torque directly. The B and 
x yy 

B gradient components are intimately related to B vv (^.B=0) . 

Z Z XX 

The By Z gradient component was exploited in Ref. 5 with the 

Spanwise Magnet roll control scheme. This scheme involves the 
addition of magnetic cores magnetized largely in the spanwise 
direction and distributed symmetrically about the model's plane 
of symmetry (Fig. 5.1). With the "2+6" electromagnet 
configuration, the required field gradient can be generated by 
the six control electromagnets, particularly effectively, in 
fact, by the central pair. 


Geometry 

Current 

5 

sense 

6 

By (T) 

B y Z ( T / m ) 

EE 

FF 

+ 

— 

0.0336 

0.0516 

0.0466 

0.0658 

Geometry 

Current 
1 2 

sense 
3 4 

By (T) 

B yz (T/m) 

EE 

FF 

+ + 

- 

0.0126 

0.0212 

0.0042 

0.0068 


5.2 Yaw angle requirement 

The difficulties inherent in the "magnetic" approach to this 
question have already been discussed. Detailed analysis would be 
a lengthy process and is quite beyond the scope of the report. 
However, modifications to the electromagnet configurations already 
studied can be anticipated. The "2+6" configuration can, it 
seems, be modified somewhat as shown in Fig. 5. 2, forming 
effectively two overlayed systems, offset at 90 degrees to each 


26 



other. It is expected that this arrangement would function 
satisfactorily around 0, 90, 180, 270 degrees yaw, but careful 
study of intermediate angles is required. 

There is a less elegant but workable alternative. A single 
point levitation system could be mounted on a large turntable, 
with the axis of rotation vertical, and simply aligned 
mechanically in the desired yaw orientation. This would easily 
give 360 degree capability. Some small magnetic maneuvering 
capability would exist around each local datum orientation. There 
are several difficulties however : 

1) Electromagnet weight. A vague estimate of the weight of 

electromagnets, based on Section 4, would be around 15 
tonnes (8 tonnes electromagnets, 7 tonnes support 
structure) . This is not an unreasonably heavy assembly 
to move with precision; astronomical telescopes are 
surely of comparable weight. 

2) Current leads and coolant lines. It is felt that careful 

design could accomodate the numerous high current 
supply leads to the electromagnets. These would be 
normally conducting where flexibility was required. 
Flexible helium transfer lines are technically feasible 
and are used in some applications, but tend to be 
expensive and delicate. 

3) The location of the electromagnet assembly would have to 

be accurately monitored. This appears straightforward. 
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6 OPERATIONAL CONSIDERATIONS 


6.1 Stray magnetic fields 

This is an important issue is the choice of site for an ML 
facility. The system creates large regions of intense magnetic 
fields. There are known biological effects of steady magnetic 
fields, though the "safe" levels of exposure are ill-defined at 
present. Several nuclear physics experiments have been operated 
on the following guidelines (Ref. 13) : 


Time Body area 


Hours Whole body 

Arms & hands 

Minutes Whole body 

Arms & hands 


Field strength (Tesla/Gauss) 


0.02 / 200 
0.2 / 2000 
0.2 / 2000 
2 / 20,000 


Several recent reports concerning biological effects of 
magnetic fields have been discovered and have been added to the 
reference list. 

Calculated magnetic field from a single optimized 
electromagnet from Section 4.1.2 is shown in Fig. 6.1. This shows 
that personnel would be excluded from the ML area during system 
operation. Since the loss of control of the levitated model is 
always a possibility (see Section 6.2) this may not be a serious 
restriction. 

More troublesome could be the effect of weaker fields on 
sensitive equipment in the vicinity of the ML system. With the 13 
inch MSBS, fields of the order of 1 Gauss have been shown to have 
an observable effect on standard CRT displays (picture 
distortion) . A 1 Gauss field criteria would only be satisfied 
around 15/20 meters from the center of the single elctromagnet of 
Fig. 6.1. 

Some form of magnetic shielding is possible, perhaps 
desirable. This would take the form of high permeability (steel, 
permalloy, mumetal) lining of an ML laboratory, such that the 
return flux paths would tend to be through the shielding rather 
than through outside air-paths. Careful calculation of the 
required shielding thickness and performance would be necessary. 


6.2 Reliability 

Extensive consideration of this point has been devoted to 
wind tunnel MSBSs. It had been considered that a Large MSBS would 
have to be very reliable, insofar as loss of control of a 
suspended model would have to be a very rare event. This is due 
to the fact that wind tunnel models tend to be expensive to 
construct, rather delicate and loss of control with wind on would 
inevitably result in the model being blown down the tunnel and 
damaged or destroyed. The difficulty in achieving high 
reliability with MSBSs arises because the inclusion of redundant 
hardware in the electromagnet array is expensive. Nevertheless, 
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it had been concluded (Ref. 5) that effective hardware redundancy 
could be achieved at reasonable cost, with careful design and 
operational techniques. 

Reliability requirements for ML may be less severe. The 
lower loads and lower dynamic force capability imply that loss of 
control might be a less violent event than is sometimes the case 
with MSBSs. Rapid shutdown (current dump) of electromagnets would 
result in the model merely falling to the floor. If the model 
were robustly made and the floor were covered with some form of 
energy absorbing material, such an event might be perfectly 
acceptable. 
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7 DISCUSSION 


Calculations performed so far seem to confirm the technical 
feasibility of the demonstration Magnetic Levitation system. 
Electromagnet and power supply specifications appear quite 
reasonable. A key design choice, inviting an early decision, has 
been identified, relating to the problem of the yaw angle 
requirement : 

a) The "magnetic" approach requires considerable design 

analysis at an early stage, but provides an elegant and 
versatile solution. 

b) The "mechanical" approach is guaranteed to work but is 

somewhat inelegant and may represent no cost saving 
over a) above. 

Whatever the choice, a comprehensive system simulation 
effort should be undertaken, in order to better identify certain 
design parameters, such as power supply voltage requirements. 

All calculated electromagnet specifications should be 
regarded as preliminary . Considerable further analysis is 
required, particularly involving the effect of magnetic couplings 
and simultaneous provision of magnetic force and moment 
components. It is certain that substantial optimization is 
possible. 
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APPENDIX 1 


The computer program FORCE (Ref. 5) was developed from the 
MIT program TABLE, specifically for the purpose of analysis of 
electromagnet geometries for MSBSs. Arrays of air-cored 
electromagnets are represented by straight line conductor 
segments and magnetic fields can then be calculated by summation 
of the effect of each segment, by repetitive use of the Biot- 
Savart law. If the magnetization of the model cores is known 
(permanent magnets) then these cores can be represented as arrays 
of magnetic dipoles. Forces and moments can now be found by 
calculating the magnetic field and field gradient components, due 
to the external electromagnets, at each dipole and summing these 
elemental forces and moments over the whole model. 

The accuracy of these methods is adequate for most MSBS/ML 
design purposes, but is dependent on the level of discretization 
of electromagnet and model cores. The serious shortcomings of 
FORCE are the lack of graphical I/O and the inability to handle 
iron electromagnet or model cores. 

A few programs designed to handle the non-linear problems of 
iron cores do exist. Notable amongst these are GFUN and TOSCA. 
GFUN relies on an integral technique and has been successfully 
used for MSBS rolling moment calculations (Ref. 5). TOSCA uses a 
differential equation method. Each program has its own advantages 
and disadvantages, but both are very sophisticated and require 
considerable care in operation, though are capable of quite 
accurate results. 


' 1. Report No. 2. Government Accession No. 

NASA CR- 178301 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Technical Background for a Demonstration Magnetic 
Levitation System 

5. Report Date 

May 1987 

6. Performing Organization Code 

7. Author(s) 

Colin P. Britcher 

8. Performing Organization Report No. 

10. Work Unit No. 

505-61-91-01 

9. Performing Organization Name and Address 
Old Dominion University 

Department of Mechanical Engineering and Mechanics 
Norfolk, VA 23508 

11. Contract or Grant No. 
NAG1-716 

13. Type of Report and Period Covered 

12. Sponsoring Agency Name and Address 
National Aeronautics and Space Administration 
Washington, DC 20546-0001 

Contractor Report 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Technical Monitor: Richmond P. Boyden, NASA Langley Research Center, Hampton, Virginia 


16. Abstract 


A preliminary technical assessment of the feasibility of a demonstration Magnetic 
Levitation system, required to support aerodynamic models with a specified clear 
air volume around them, is presented. Preliminary calculations of required sizes of 
electromagnets and power supplies are made, indicating that the system is practical. 
Other aspects, including model position sensing and controller design, are briefly 
addressed. 


17. Key Words {Suggested by Author(s)) 


18. Distribution Statement 


Magnetic Levitation 
Magnetic Suspension 


Unclassified - Unlimited 
Subject Category - 09 

19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price 

Unclassified 

Unclassified 


37 

A0 3 


For sale by the National Technical Information Service, Springfield, Virginia 22161 

















